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Abstract 



The interaction energy between two black holes at large separation 
distance is calculated. The first term in the expansion corresponds to the 
Newtonian interaction between the masses. The second term corresponds 
to the spin-spin interaction. The calculation is based on the interaction 
, energy defined on the two black holes initial data. No test particle approx- 

• imation is used. The relation between this formula and cosmic censorship 

fSJ ' is discussed. 

o ■ 

qhI 1 Introduction 



The purpose of this article is to prove, under appropriate assumptions, the fol- 
lowing statement: the interaction energy, at large separation distance I, between 
two black holes of masses Mi , M2 and spins Ji , J2 is given by 



. -M1M2 , -Ji- J2 + 3(Ji-n)(J2-n) ... , ^ 

H , E = ; 1 — — + higher order terms, (1) 

I l'^ 

where n is a unit vector which points inward along the line connecting the black 
holes. Before giving a precise definition of the parameters involved in Eq. (|l|), 
I want to discuss its physical meaning. 

The first term in Eq. (|^) has the Newtonian form. For two point particles 
of masses Mi and M2 separated by an Euclidean distance the Newtonian 
interaction energy between them is given by —MiM2/l- The fact that this 
term appears also for a two black holes system in General Relativity can be 
expected from the weak field limit of Einstein's equations. The second term 
in Eq. (||), which involves the spins, has analogous form to the dipole-dipole 
electromagnetic interaction; there exists an analogy between magnetic dipole in 
electromagnetism and spin in general relativity (see ]35[|). In the electromagnetic 
case, / is the Euclidean distance between two charge distributions and Ji , J2 the 
corresponding dipole moments of them. However, the gravitational black hole 
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spin-spin interaction has the opposite sign to the electromagnetic one. The first 
evidence of this fact was given by Hawking . I want to reproduce Hawking's 
argument here because it points out the connection between Eq. (|I]) and the 
cosmic censorship conjecture (see also the discussion in Q). In the argument, 
we assume the following two consequences of weak cosmic censorship and the 
theory of black holes (cf. ^ |37), see also |3§1): 

(i) Every apparent horizon must be entirely contained within the black hole 
event horizon. 

(ii) If matter satisfies the null energy condition (i.e. if Tabk^kf' > for aU nuU 
k"'), then the area of the event horizon of a black hole cannot decrease in 
time. 

We also assume: 

(iii) All black holes eventually settle down to a final Kerr black hole. 

Consider a system of two black holes such that, at a given time, the separa- 
tion distance between them is large. Then, there must exist a Cauchy surface in 
the asymptotically flat region of the space time such that the intersection of the 
hypersurface with the event horizon has two disconnected component of areas 
Ai and A2. Since the black holes are far apart, these areas can be approximated 
by the Kerr formula 

Ai = 8n (^M^ + ^JMf-J^^ , ^2 - Stt (^A/| + y^Af| - . (2) 

At late times, after the collision, the system will settle down to a Kerr black hole. 
Hence, there must exist another Cauchy hypersurface such that its intersection 
with the event horizon will have area 

Af = 8TT (AI^j + ^Mj- jfj , (3) 

where Mf is mass of the final black hole and Jf is its final angular momentum. 
By (ii) we have 

Af>Al^-A2. (4) 

Since gravitational waves have positive mass, we also have 

Mf<Mi+M2. (5) 

In general, gravitational waves will carry angular momentum. But in axially 
symmetric space-times the total angular momentum is a conserved quantity, 
since it can be defined by a Komar integral (cf. ||2^ and also |^). Then, in 
this case we have 

Jf = Ji + J2. (6) 
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Using Eqs. (^), (^), (^) and (||) it is possible to obtain an upper bound, which 
depends on Ji and J2, to the total amount of radiation emitted by the system 
Ml + M2 — Mf. It can be seen that if Ji and J2 have the same sign, this 
upper bound is smaller than if they have opposite sign. This suggests that 
there may be a spin-spin force between the black holes that is attractive if the 
angular momentum have opposite directions and repulsive if they have the same 
direction. Presumably, in the second case the system expends energy in doing 
work against the spin repulsive force, and for this reason this energy is not 
available to be radiated via gravitational radiation. 

Hawking's argument only suggests that the spin interaction energy between 
black holes has in fact this sign dependence with respect to the spins. It is not a 
proof, first because there is no proof for the weak cosmic censorship conjecture 
(i)-(ii) and for the assumption (iii). Second, because even if we assume (i)- 
(iii) the argument only shows that an upper bound of the total amount of 
radiated energy has this sign dependence in terms of Ji and J2, but the real 
amount of gravitational radiation can, in principle, have other dependence. In 
fact, the total amount of gravitational radiation produced by such systems, as 
numerical studies show, is much smaller than this bound. This upper bound 
is 50% of the total mass when the spins are antiparallel, the black holes are 
extreme (J^ = M), and have equal masses; when the spins are zero or when the 
black holes are extreme with parallel spins, the upper bound is 29% of the total 
mass. On the other hand, in the numerical calculations the maximum amount 
of radiation emitted by this type of system is about 3% of the total mass, see 

j|] for a recent calculation and also for an up to date review on the 
subject. However, the numerical studies show that the system indeed radiates 
less when the spins are parallel than when they are antiparallel. Moreover, Wald 
[ p5| proves that the interaction energy between a test particle with spin J2 and 
a stationary background of spin Ji has precisely this sign dependence. Wald 
shows that the spin-spin interaction energy has the form 

-Ji ■ J2 + 3( Ji • n) (Ji ■ h) 



where I and n are defined as follows. The stationary field is expanded at large 
distance with respect to Cartesian asymptotic coordinates a;*, here / is the Eu- 
clidean radius with respect to and = /I. Eq. (0) has been also proved 
by D'Eath using post-Newtonian expansions jl^. It is important to note that 
Eq. (0) gives an indirect evidence in support of (i)-(iii). 

In this article I want to prove Eq. (0) without using cither the particle or 
post-newtonian approximation. The proof is based on an interaction energy 
defined on the two black hole initial data. This interaction energy is genuinely 
non linear; it does not involve any approximation. 

The plan of the paper is as follows. In section the main results are given. 



In section theorem 2.2 is proved; in section M we prove corollary 2.3. Finally, 



in section ra an alternative definition of the interaction energy is discussed. 
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2 Main Result 



The strategy I will follow was given by Brill and Lindquist . It is based on 
the analysis of initial data set with many asymptotic ends. An initial data set 
for the Einstein vacuum equations is given by a triple (S, hab, Kab) where S* is a 
connected 3-dimensional manifold, hab a (positive definite) Riemannian metric, 
and Kab a symmetric tensor field on S. They satisfy the vacuum constraint 
equations 

^kab - Dak = 0, (8) 

R + k^- kabk'^' - 0, (9) 

on S, where Da is the covariant derivative with respect to hab, R is the trace 
of the corresponding Ricci tensor, k = h°'^kab, and a,b,c... denote abstract 
indices. Tensor indices of quantities with tilde will be moved with the metric 
hab and its inverse h""^. The data will be called asymptotically flat with A'^ + 1 
asymptotic ends, if for some compact set 51 we have that S\Q — X]^o ^k, where 
Sk are open sets such that each Sk can be mapped by a coordinate system 
diffeomorphically onto the complement of a closed ball in R^^ such that we have 
in these coordinates 

kj^il + ^)S..,+0{r-'), (10) 

k,, = 0(r-2), (11) 

as r = (X^j^iC^^"*)^)^^^ ^ oo in each set Sk; where i,j---, which take values 
1,2,3, denote coordinates indices with respect to the given coordinate system 
, and Sij denotes the flat metric. We will call the coordinate system a;* an 
asymptotic coordinate system at the end k. Each asymptotic region Sk has a 
different asymptotic coordinate system. The constant denotes the ADM 
mass[|l| of the data at the end k. These conditions guarantee that the mass, the 
linear momentum, and the angular momentum of the initial data set are well 
defined at every end. 

For A^ > 1, this class of data contains, in general, apparent horizons. The 
existence of apparent horizons leads us to interpret these data as representing 
initial data for black-holes. Their evolution will presumably contain an event 
horizon, according to the standard theory of black holes ||2^. The validity of 
this picture depends, of course, on the cosmic censorship conjecture. The only 
statement about the evolution of the data that we can make is the geodesic 
incompleteness of the space time. In general, in order to prove the geodesic 
incompleteness of a space time, one needs to know that the data contain a 
trapped surface in order to apply the singularities theorems |2^. However, in 
this particular case, since the topology of the data is not trivial, the geodesic 
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50 , Mo , Jq 



Figure 1: Initial data with three asymptotic ends (TV = 2). For each asymptotic 
region Sk, we have the corresponding mass Mk and total angular momentum 



incompleteness of the space time follows directly from a theorem proved by 
Gannon 

For simplicity we will fix iV = 2, see Fig. ^ In this case the data can be 
interpreted as initial data with two black holes. This interpretation is suggested 
by the following fact: when an appropriate distance parameter is large compared 
with the masses Mk, then it can be seen numerically that only two disconnected 
apparent horizons appear. For time symmetric data, these numeric calculations 
have been done in ; the non-time symmetric case has been studied by Cook 
(see [|4| and references therein) . It is not clear that the number of apparent hori- 
zons is the number of black holes contained in the data, since even when there 
are two disconnected apparent horizons, the intersection of the event horizon 
with the initial data can be connected. However, at large separation distance, 
this seems to be a reasonable assumption, which is confirmed by the numerical 
evolutions [ p9| . 

Brill and Lindquist define the following interaction energy at the end k 



The energy Ek is a geometric quantity; its definition does not involves any ap- 
proximation. The question now is how to calculate Ek in terms of physically 
relevant parameters. The first problem is how to define an appropriate separa- 
tion distance between the black holes. When there are two apparent horizons, 
there is a well defined separation distance Ij^ defined as the minimum geodesic 
distance between any two points in the two different horizons, see Fig. |^. 

However, the distance Ij^ is hard to compute. The location of the apparent 
horizons can be calculated only numerically. Since we are only interested in 
the energy at large separations, instead of If^ we will use another parameter I, 
and we will argue that « Z in this limit. The definition of the parameter 
I is related to the way in which one can construct solutions of the constraint 
equations with many asymptotic ends. The conformal method (cf. pl|, p2| 




N 




(12) 



k'=0 
k'^k 
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Apparent 




Figure 2: An initial data set with three asymptotic end points and only two 
disconnected apparent horizons of area Ai, A2, and radii R2. The points 
ii and 12 represent the two other infinities 1 and 2. The geodesic distance is 
computed with the physical metric hab- The parameter I is computed with the 
conformal metric. The geodesic distance between ii and 12 with respect to the 
physical metric is infinite. 



and the references therein) is a general method for constructing solutions of 
the constraint equations. We assume that hab is a positive definite metric with 
covariant derivative D^, and K'^^ is a trace- free (with respect to hat), symmetric 
tensor, satisfying 

Da.K''^ = Q on S. (13) 

Let y be a solution of 

ih<^ - -^i^afcif'V"^ on S, (14) 

o 

where Lh = D'^Da — R/8 and R is the scalar curvature of the metric hat- Then 
the physical fields {h,K) defined by Kb = <P^hab and 7^"'' = Lp-^^K'''' wiU 
satisfy the vacuum constraint equations on S. We have assumed that K"''^ is 
trace- free; hence K"-^ will be also trace- free with respect to hab- That is, the 
initial data set will be maximal. 

To ensure asymptotic flatness of the data at the each end fc, we will require 
the following boundary conditions. Let i\ and i2 be two, arbitrary points in M^, 
with coordinates x\ and in some Cartesian coordinate system a;*. Define the 
manifold 5 by 5 = \ {ii, Z2}. Assume that hab is regular on M^. At infinity 
we will impose the following fall off behavior 

hij=5ij+0{r-^), (15) 
K"^ = 0{r-^), (16) 
<f = l + 0{r-^). (17) 
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At the points ii and «2 we require 

K^'^Oir^^), K^'^0{r^% (18) 



where 




and 



3 \ 1/2 



Um riLp — — lim r2'^ = — ^, (20) 

r-i^O 2 ' r2— »0 2 



where mi and m2 are positive constants. Note that both (f and K""^ are singular 
at ii, 12- 

One can prove that the data so constructed will be asymptotically flat at the 
three ends. We have made an artificial distinction between the end 0, given by 
r ^ cxD, and the ends 1 and 2. It is possible to discuss the same construction in 
a more geometrical way, such that all ends are treated equally; see [|], pl| , 
[ p^ . However, since our final goal is to calculate the interaction energy at one 
end, it is convenient to make this distinction. The coordinate system and 



the corresponding flat metric in the expansion Eq. (15), gives the Euclidean 
distance I between ii and 12 

/ 3 \ 1/2 

/=(E(^2-^in , (21) 

which will be our separation distance parameter, see Fig. |^. 

In general, Eq. (|lj) is non-linear. However if we assume that the data is 
time symmetric, i.e. if"'' = 0, then it becomes a linear equation for Lp. If we 
assume that the conformal metric is flat, we obtain a Laplace equation for if. 
The solution of this equation that satisfies the boundary conditions Eqs. ( pO| ) 
and ( [TtI ) is given by 

This solution was found by Brill and Lindquist in . In this case it is possible 
to calculate explicitly the interaction energy ( p^ in terms of the masses and 
the separation distance. The result is the following. 

Theorem 2.1 (Brill-Lindquist) Let hat the flat metric and K""^ — 0. Then 
the interaction energy defined by Eq. is always negative. Moreover, when 

I is large compared with Mk the following expansion holds 

Eq = ~7~~~ + higher order terms. (23) 
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Giulini has computed the higher order terms for these data and other con- 
formally flat time symmetric data with different topologies. In those examples 
the Newtonian term is invariant but the higher order terms depend on the par- 
ticular initial data. 

In order to discuss spin-spin interaction, we need initial data with non trivial 
angular momentum, that is we have to allow for non trivial extrinsic curvature 
in the data. At each end we have the angular momentum Jk given by 

Jl = ^Yim f rKbcv}'t^'"'nddSr, (24) 

where Sr is a two sphere defined in the asymptotic region Sk and n° is its 
outward unit normal vector. In Eq. ( p^ we can use either K""^ or K""^ because 
the conformal factor satisfies In general the angular momentum at each 

end is not determined by the intrinsic angular momentum of each black hole. It 
includes also the angular momentum of the gravitational field surrounding the 
black holes. Then, in general, there is no relations between Jq, Ji and J2, these 
three quantities can be freely prescribed. But in the presence of symmetries 
these quantities can not be given freely any more. Moreover, in the presence 
of conformal symmetries of the metric there exists a well defined quasilocal 
definition of angular momentum. Assume that is a conformal Killing vector; 
that is, a solution of the equation {ChS)""^ = 0, where 

[ChiT'' = D'^C'' + D^C - I h"^ D^^". (25) 

If the initial data is maximal, i.e., K — 0, then the vector K'^^'^i, is divergence 
free. Hence, for each conformal symmetry we have the associated integral 

1^= f K^^bTiadS, (26) 

JS 

where S* is a close 2-surface and n" its outward unit normal vector. This integral 
is a conformal invariant. It can be calculated also in terms of tilde quantities. 
The integral (p6| ) will be non zero only if the vector K'^'^^i, is singular at some 
points; in our case it will be singular at two points: the location of the holes. 
Then the integral in Eq. ( p6| ) will have three different values depending on 
whether the surface S encloses one hole, two holes or no hole. In the later case, 
=0. If we chose to be a rotation, will gives the corresponding component 
of the quasilocal angular momentum. If the data is conformally flat we have 10 
conformal Killing vectors. In particular, we have the three rotations and hence 
the complete definition of quasilocal angular momentum. These quantities will 
be defined only on this slice and will generally not be preserved in the evolution. 
They will be only preserved if the space time admits a Killing vector. In this 
case they will coincide with the corresponding Komar integral. The space time 
will admit a Killing vector field if i^" is a Killing vector for the whole initial 
data; that is, £^hab ~ £^K°''' = 0, where £^ is the Lie derivative with respect 
to A conformally flat, maximal, slice can be interpreted as a instant of time 



8 



in which the gravitational field carries no angular momentum and no linear 
momentum itself, and hence these quantities are carried only by the "sources" , 
which in this case are the black holes. Data containing matter with compact 
support can be also constructed. 

There exist in the literature other definitions of quasilocal angular momen- 
tum Q @@), which are apphcable for an arbitrary closed 2-surfaces in 
the spacetime. It is not clear if any of these definitions will agree with Eq. ( p6[ ) 
in the particular case of 2-surface lying on a conformally flat 3-hypersurface. 

From the discussion above, we conclude that in the case of conformally flat, 
maximal data we have, 

Jo° + J1 + - 0. (27) 

For an observer placed in the asymptotic end the system will look like two 
black hole with spins — J° and — J2 , and the total angular momentum will be 
Jq = —Jf — J£. For a more general discussion of conformal symmetries on 
initial data see [pj and JTzt ; in particular in those articles a generalization of Eq. 
( p7| ) which includes linear momentum is proved. 

Bowen and York obtain a simple model for a conformally flat data set which 
represents two black holes with spins . Brandt and Briigmann study these 
data with the + 1 asymptotic ends boundary conditions given by Eqs. (p!8[), 
I6h, (EOh and (O). For these data, the conformal second fundamental form is 



j^ab ^ j^ab j^ab^ ^28) 



T^ab _ 6 „(<»,b)cd j „ j^ab _ ^ J>)cd j „ (r,Q.\ 



given by 



where 



and 



(30) 



ri r2 

where Jf , J| are constants and e^^** is the flat volume element. One can check 
that the constants Jf and J| give the the angular momentum at the ends 1 and 
2 respectively. The tensors ( p9| ) are divergence free and trace free with respect 
to the flat metric in M.^ \ {ii,i2}- 

The first result of this article is the following theorem. 

Theorem 2.2 Let hab be the flat metric and K"-'' be given by Eqn. (p^j. Then, 
the interaction energy ^T^ ) is given by 

Eq ~ ^ H ^— !- — — l2L_t — 11 -|- higher order terms. (31) 
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The expansion of the dimensionless quantity Eo/l is made in terms of the di- 
mensionless parameters Mi/Z, M2/I, Ji/P and J^jl'^- By higher order terms 
we mean terms of cubic order in those parameters. We prove Theorem [j.2| in 
section ^. 

Using similar arguments, Gibbons ||2^ has obtained the quahtative sign 
dependence of the spin-spin interaction for a certain class of axially symmetric, 



conformally flat, data. Note that in theorem 2.2, the data are non axially sym- 
metric in general, since the spins can point in arbitrary directions. Bonnor [Q 
has studied the spin-spin interaction using an exact, axially symmetric solution 
of the Einstein-Maxwell equations, his result also qualitatively agrees with the 



spin-spin term in Eq. (31) 



The Bowen-York data are, of course, very special. The natural question is 
how to generalize theorem ^.2| for more general data. For general asymptotically 
flat data with three asymptotic ends we can not even expect to recover the 
Newtonian interaction term. Take for example a time symmetric initial data 
with only two ends. Choosing the conformal metric appropriately, one can 
easily construct data such that the difference M\ — Mq is arbitrary. That means 
that we are putting more radiation in one end than in the other. If we add a 
third end with small mass, the new interaction energy will be dominated by the 
difference Mi — Mq and hence will be not related with any Newtonian force. 



Hence, theorem 2.1 is not true for general asymptotically flat metrics with many 
asymptotic ends. 

The interaction energy defined by (|l^) can have the meaning of a two black 
holes interaction energy only if is it possible to distinguish in the data two objects 
that are similar to the Kerr black hole when the separation distance is large. I 
will call this class of data two Kerr-like black hole initial data. The existence 
of these data has been proved in . The following two properties of the 

Kerr initial data, in the standard Boyer-Lindquist coordinates, are important: 
i) The data are conformally flat up to order O(J^) ii) The leadingorder term of 
the second fundamental form is given by the Bowen-York one (^9|). Then one 
can expect that Eq. (|3^) is unchanged in the principal terms for this class of 
data. 

More precisely, a two Kerr-like black hole initial data can be constructed 
as follows (see [|l^ and ||l^ for details). Take a slice of the Kerr metric, with 
parameters Mk^ and Ji, in the standard Boyer-Lindquist coordinates. Choosing 
the appropriate conformal factor, the conformal metric can be written in the 
following form 

hab=3at + h!^^, (32) 

where h^^ ^ 0{Jf). In the same way the conformal second fundamental form 
can be written as 

K'^\ ^ Kf + K'^\, (33) 
where Kfj'^ = 0(Ji ) and Kf^ is given by (p9|). Take another Kerr metric, with 
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parameters and J2, and define the following conformal metric 

h^'' = 5a, + h"^t+h^S^ (34) 
and the following conformal second fundamental form 

K-^K = Kt + Kt + i^h^.^r", (35) 

where the bar means the trace free part of the tensor with respect to the metric 
(HI) and is chosen such that Kkk divergence free and trace free with 
respect to the metric (^). I n [l6| p5[ it has been proved that such vector w°' 
exists and is unique. Using (^4[) and (^5|), solve Eq. ( p^ with the boundary 
conditions (|2^) and (^^ where 

m,^^M]~Jl/Ml^, m^ = ^Ml^-Jl/Ml^. (36) 

The existence of a unique solution has been proved in Q . 

If we chose and to point in the same direction, the data will be axially 
symmetric. In this particular case, we can use the integral (|2^) to calculate the 
quasilocal angular momentum of each of the black holes. The result will be Jf 
and J2 . However, in general, this class of data will admit no conformal Killing 
vector. In this general situation, it is very hard to compute the quasilocal spins 
of each of the black holes. However, when the separation distance is large, 
and J2 will give approximately the angular momentum of each of the black 
holes because this class of data have a far limit to the Kerr initial data. In other 
words, and give the spins of one black hole when the parameters of the 
other are set to zero. 

For this class of data we have the following result. 

Corollary 2.3 For the two Kerr-like data defined above, the formula ^^\ ) for 
the interaction energy holds. 

We prove this Corollary in section ^. 



3 Interaction Energy for the spinning Bowen- 
York initial data 



In this section we will prove theorem 2.2. For Bowen-York data, with the 
boundary condition ( |20| ) and (p^, Eq. (|14|) for the conformal factor can be 
written in the following forni[^ 

= g^7 ' </3 = </3o + u, (37) 

with the boundary condition 

lim u = 0, (38) 
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where ifo is defined in Eq. (^2|), with mi and m2 arbitrary positives constants, 
and if"'' = Kf + Kf is given by (|2|). 

The coordinates are asymptotic coordinates for the end 0. The total mass 
at is given by 

Mo = mi + m2 + 2uoo, (39) 

where Uoo is the term which goes hke 1/r in the solution u of Eq. and is 
given by 

uoo-^ f ^^^dx\ (40) 
47r 7r3 8ip^ 

We want to calculate the masses Mi and M2 for the other ends. The asymptotic 
coordinates for the other ends are 



-t _ '"-1 V-^ ~ -^1) -i _ '"-2 V-^ ~ -^21 fy,■,^ 



I2 ' ^42 ~ ~1 

Take for example the end point ii, in coordinates we have 



where we have used 

r2 =; + 0(l/fi), (43) 

and f 1 is the Euclidean radius with respect to x . Then the mass at this end 
is given by 

Ml =mi (l + ^ + u(ii)) , (44) 

where u(ii) denote the value of the function u at the point ii. In analogous way 
we obtain the mass at the end 2 

/ m<i \ 
M2 = ma 1 + — !- + u{i2) . (45) 



21 

Then the interaction energy at the end iq is given by 

Eo — Mq ~ Ml ~ M2 — J h 2moo - miu(ii) - m2u{i2). (46) 

Using Eq. (^7|) and the Green's function for the Laplacian, we obtain the 
following integral representation of the terms involving u in Eq. ( p6| ) 

2woo - miu(ii) - m2u{i2) = —— / = — 1 " 7; 7;— ] ■ (47) 

167rjR3 if' \ 2ri 2r2 ) 
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The formula involves the unknown function u. Using the fact that u is 
0{J^), we make an expansion of this integral in terms of the parameters Ji/l"^, 
Jijl"^ 1 mi /I and m^jl- We obtain that the first non trivial term is given by 

luoo - miu{i^ ~ m.2u{i-2) E'' , (48) 

where 

E' = ^ f KfK2abdx. (49) 

The interaction energy, up to this order, is given by 

E = -^+E^ (50) 

where we have used Eqs. (^J) and (^ ) to replace nik by Mk, since up to this 
order they are equal. 

All that remains is to compute the integral E''. This integral can, in princi- 
ple, be calculated explictly from the expressions (|2^). However such a calcula- 
tion is very complicated. Instead of this we will calculate (|4^) in the following 
way. The tensors ( |29| ) can be written like 

{Csvif = , {Csv2r'' = Kf (51) 

where is the conformal Killing operator defined in Eq. (|2^) for the flat 
metric, and 

v\ = -^^\h,nikri\ (52) 

vl = -^^^J2,n2ur2\ (53) 

Let and B^^ be small balls centered at «i and 12 respectively, of radii ei and 
£2. We have that 

E'^^ lim / KfK2abdx. (54) 

»7r ci,c2^o Jk3-b,-^ 

Using the Gauss theorem in M.^ — B^^ — B^^ we obtain 

KfK2 abdx = -2 [ Kfvi bui a dS,, - 2 / Kfvi „ dS,, , 

(55) 

where n° and are the outward normals to the two surfaces B^^ and . In 
the limit ei ^ the first integral vanishes. We use here that Xj'' is regular 
in B^^. The second integral, in the limit 62 ^ can be easily calculated. We 
obtain 

^ -Ji ■ •^2 + 3(Ji ■ n){Ji -fi) ^gg^ 
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4 Interaction energy for the two Kerr-like black 
holes initial data 



The two Kerr-like black hole initial data set are solutions of the following equa- 
tion 

LhKKipKK = on S, (57) 

with the boundary conditions (20) and (p7|), where /i^^ and K^^j^ are given by 
(p4[), ( p^ ) and mi, m2 are given by ( |36| ) in term of the Kerr parameters. 

The conformal factor for the Kerr initial data, with parameters Mk^ , Ji can 
be written in the following form 

VKi=l + ^ + ^R,, (58) 

where (pu-^ — 0(Ji). We can decompose the two Kerr solution in the following 
form 

ifKK =VQ + VRi + + (59) 

Using that ( p8| ) is a solution for one Kerr initial data, we obtain that the first 
term in the expansion in Ji and J2 of the function u satisfies the following linear 
equation 

An « (60) 

Hence, the spin-spin interaction term has the same form as the Bowen-York 
one. In an analogous way to the previous section, we obtain that the interaction 
energy is given by 

e„.-MiMi^e: (61) 

where is given by Eq. (E^). 



5 Discussion 

We have shown, using the interaction energy defined by Eq. (|2|), that the 
spin spin interaction between black holes of arbitrary masses and spins has an 
expansion of the form (|3^). This formula has been previously derived using 
a test particle approximation and post-newtonian expansions The 
main improvement of the present calculation, beside its simplicity, is that no 
approximation is used in the definition of the interaction energy. Moreover, the 
very definition of the interaction energy involves black holes, in contrast with 
previous calculations where the black holes appear indirectly. 
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The interaction energy defined by Eq. (|T^) uses the fact that we are choosing 
a particular topology for the initial data, but other topologies are possible. 
Examples of different kind of topologies where we can not use Eq. (12) are the 
Misner topology of two isometric sheets [Q, the Misner wormhole [31| or even 
initial data with trivial topologies which contain an apparent horizon | J . If the 
initial data have k disconnected apparent horizons of area Ak , we can define the 
individual masses as follows 



Mh.=\I^- (62) 



Then, the interaction energy is given by the formula 

N 

Eh=M-Y,Mhi. (63) 

fc'=0 
k'^k 

What is the relation between Ek defined by Eq. (^2|) and Eh defined by (|63|)? 
Note that Eh is much harder to compute then Ek- I want to argue that Ek 
will presumably give the sa me r esult as Ek for the leading orders terms in the 
expansion given by theorem |2.2| . Assume that the data is such that when the 
separation distance parameter is large, then the areas can by approximated by 
the Schwarzschild formula Ak ~ 167rM| plus terms of order J^. Then the radius 
of the horizons will be Rk « 2Mfe. The distance If^ will differ from I in terms 
order 0{M,J^). Then we can replace I by in theorem and Mh^. by Mk in 
Eq. (|l2|), up to this order. 

It is interesting to note that the interaction energy Eh has been used in a 
different context, namely to determine the last stable circular orbit in a black 
hole collision [[131. 
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